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The problem of calculating collective density fluctuations in quantum liquids is revisited. A fully 
quantum mechanical self-consistent treatment based on a quantum mode-coupling theory [E. Rabani 
and D.R. Reichman, J. Chem. Phys, 116, 6271 (2002)] is presented. The theory is compared with 
the maximum entropy analytic continuation approach and with available experimental results. The 
quantum mode-coupling theory provides semi-quantitative results for both short and long time 
dynamics. The proper description of long time phenomena is important in future study of problems 
related to the physics of glassy quantum systems, and to the study of collective fluctuations in Bose 
fluids. 



I. INTRODUCTION 

The study of the dynamical properties of quantum liq- 
uids has a storied history^i The interplay of nuclear dy- 
namics, particle statistics, dimensionality, disorder and 
temperature can lead to, or suppress, dramatic effects 
such as superfluidity in "^He and superconductivity in the 
superfluid Fermi liquid state of conducting electrons in 
metalsi^ Due to the great computational difficulties pre- 
sented by a direct assault on the real-time many-body 
Schroedinger equation, microscopic approaches to cal- 
culate the frequency dependent response in condensed 
phase disordered quantum systems are approximate and 
generally rely on somewhat uncontrolled approximations. 

One approach that has been useful in a variety of 
physical contexts is the analytic continuation of numer- 
ically exact imaginary time path-integral data."^'^ This 
approach has been fruitful for the study of quantum im- 
purity models such as the Anderson Hamiltonian,^ as well 
as models of correlated electrons such as the Hubbard 
Hamiltonian.^ Berne and coworkers have successfully ap- 
plied the maximum entropy (MaxEnt) version of analytic 
continuation to study the dynamics of electrons in simple 
liquids,^'® vibrational and electronic relaxation of impu- 
rities in liquid and solid hosts, ^'^^ and adiabatic rate con- 
stants in condensed phase environments^ii^ Boninsegni 
and Ceperley have applied MaxEnt to study the dynamic 
structure factor of liquid ''He above and below the A- 
transition.-'^'^ For the normal state, they find good agree- 
ment with experimental results, while the agreement is 
significantly worse in the superfluid state. In particular, 
sharp quasiparticle peaks are not well resolved with the 
MaxEnt approach. ^■^ This failure arises from the intrin- 
sic ill-posed nature of the numerical continuation to the 
real-time axis. The conditioning of the data which helps 
alleviate numerical instabilities makes the differentiation 
of fine energy scales difficult, resulting in smoothing of 
the response functions. 



In several recent papersraiiS*iSiiLi& we have explored a 
molecular hydrodynamic approach for calculating the dy- 
namic response functions of quantum liquids. In partic- 
ular, we have formulated a detailed generalization of the 
classical mode-coupling approach to the quantum case^ 
We have studied both single particle as well as collec- 
tive fluctuations, with good agreement with recent ex- 
periments performed on liquid para-hydrogeni2iS2i2iiS2iS 
and liquid ortho-AewteTminMi^^^ Unlike the MaxEnt 
technique, our approach is a theory and thus provides 
additional insight into the dynamics of quantum liquids. 
This fact implies several additional advantages. First, be- 
cause our theory is the direct analog of the approach used 
in the study of quantum spin glasses ^L^SiSS interesting 
connections between quantum systems with and without 
quenched disorder may be made. For example, our recent 
study of the spectrum of density fluctuations in liquid 
para-hydrogen, where finite frequency quasiparticle peaks 
appear in the dynamic structure factor, S{q, lj), resonates 
with the general result of Cugliandolo and Lozano that 
"trivial" quantum fluctuations add coherence to the de- 
cay of dynamics correlators at short-times.^Q Further- 
more, in principle we could study a range of problems 
related to the physics of glassy quantum systems, includ- 
ing the aging behavior of an out-of-equilibrium quantum 
liquid. Because the MaxEnt approach relies im- 
plicitly on the equilibrium formulation of quantum sta- 
tistical mechanics, this approach cannot be used to in- 
vestigate such questions. 

In order to study such interesting problems, as well as 
the challenges imposed by nontrivial collective behavior 
exhibited, for example, by superfluid ^He, we need to 
further develop our theoretical apparatus. In our pre- 
vious studies, we imposed self-consistency in the study 
of single-particle properties, but not in the more dif- 
ficult case of collective properties such as the dynam- 
ics of density fluctuations In this work, we make 
several accurate and physically motivated approxima- 
tions that render the formal equations presented in Ref. 
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amenable to direct numerical investigation. We re- 
visit the problem of the study of S(q, u>) in liquid ortho- 
deuterium and para^hydrogen, and show that a fully self- 
consistent treatment of density fluctuations greatly im- 
proves the agreement between the low frequency behavior 
seen in experiments and corrects the behavior of our pre- 
viously used quantum viscoelastic theoryji^iiii^ This low 
frequency behavior associated with the long time relax- 
ation of density fluctuations in these liquids is not well 
captured by the MaxEnt approach. 

Our paper is organized as follows: In section^]we pro- 
vide an overview of our self-consistent quantum mode- 
coupling approach to density fluctuations in quantum 
liquids. Furthermore, we discuss the improvements of 
the present approach and the physical approximation 
that are introduced to make our current study amenable 
to path integral Monte Carlo (PIMC) simulation tech- 
niques. In section ITTll we discuss the MaxEnt approach 
used in the present study. Results for collective den- 
sity fluctuations in liquids orf/io-deuterium and para- 
hydrogen is presented in section IIVI We compare the 
predictions of our quantum mode-coupling theory to the 
MaxEnt results and also to available experimental re- 
sults. Finally, in section El we conclude. 



II. SELF-CONSISTENT QUANTUM 
MODE-COUPLING THEORY 



We begin with the definition of two dynamical vari- 
ables, the quantum collective density operator 



N 



(3) 



and the longitudinal current operator 



2m\q\ 



N 

E 



[(q-Pa)e^'»'°+e^''^°(p..q)], (4) 



where f „ is the position vector operator of particle a with 
a conjugate momentum pa and mass m, and N is the 
total number of liquid particles. The quantum collective 
density operator and the longitudinal current operator 
satisfy the continuity equation = iqjq, where the dot 

denotes a time derivative, i.e. i^— jr[H,pq]. 

Following the projection operator procedure outlined 
in our early work,^® the time evolution of the Kubo trans- 
form of the intermediate scattering function, F'^{q,t) — 
jj-{p\, Pq{t)) ((•■■) denote a quantum mechanical ensem- 
ble average), is given by the exact quantum generalized 
Langevin equation 

F^iq,t)+u;liq)F^iq,t) 

ft . (5) 

+ / dt'K''{q,t~t')F''{q,t') = 0, 



In this section we provide a short overview of our quan- 
tum mode-coupling approach suitable for the study of 
collective density fluctuations in quantum liquids. For 
a complete derivation of the equations described below 
the reader is referred to Ref. [l^- The formulation of 
the quantum mode-coupling theory (QMCT) described 
in Ref. [3 is based on the Kubo transform of the dynami- 
cal correlation function of interest. For collective density 
fluctuations the experimental measured quantity is the 
dynamic structure factor, S{q, lo), which is related to the 
Kubo transform of the intermediate scattering function, 
F-{q,t) by: 



S{q,uj) 



coth 



+ l\S-{q,Lo) (1) 



where the Kubo transform of the dynamic structure fac- 
tor is given in terms of a Fourier transform of the Kubo 
transform of the intermediate scattering function: 



dte''^^F''[q,t). 



(2) 



In the above equations (3 — is the inverse tempera- 
ture, and the superscript k is a shorthand notation of the 
Kubo transform (to be described below). The quantum 
mode-coupling theory described in this section provides a 
set of closed, self-consistent equations for F'^{q,t), which 
can be used to generate the dynamic structure factor us- 
ing the above transformations. 



where the notation k implies that the quantity under 
consideration involves the Kubo transform given by^i 



1 



(6) 



and H is the Hamiltonian operator of the system. In 
the above equation a'^((?) is the Kubo transform of the 
frequency factor given by: 



1 



(7) 



with S'^{q) = F^{q,Q) = ^(pJj,/5^(0)) being the Kubo 
transform of the static structure factor. The Kubo 
transform of the memory kernel appearing in Eq. (O 
is related to the Kubo transform of the random force 



Rn, 



dt 



SHq) 



K-{q,t) 



/5q, and is formally given by 



(i?t,e^(i-^")^*i?«). 



NJ'^{q) 



(8) 



where J'^{q) — j^{jq,jqiO)) is the Kubo transform of 
the zero time longitudinal current correlation function, 
C — [H, ] is the quantum Liouville operator, and 

'P'^ — ■^^ta^'A" is the projection operator used to 
derive Eq. jSJl with the row vector operator given by 
A = (Pq:iq) [sGc Ref. for a complete derivation of 
Eqs. ®, O, and ©]. 
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The above equations for the time evolution of the 
intermediate scattering function is simply another way 
for rephrasing the quantum Wigner-Liouville equation. 
The difficulty of numerically solving the Wigner-Liouville 
equation for a many-body system is shifted to the difh- 
culty of evaluating the memory kernel given by Eq. JHJ. 
To reduce the complexity of solving Eq. © we follow the 
lines of classical mode-coupling theory'^^-^'^''^'* to obtain 
a closed expression for the memory kernel. A detailed 
derivation is provided elsewhere^ 

We use the approximate form K'^{q,t) — K'J{q,t) + 
K!^{q,t) where the "quantum binary" portion K'^{q,t) 
and the quantum mode-coupling portion K'^{q, t) are ob- 
tained using the standard classical procedure r^i^^iSi but 
with a proper quantum mechanical treatment. The ba- 
sic idea behind our approach is that the binary portion 
of the memory kernel contains all the short time infor- 
mation and is obtained from a short time moment ex- 
pansion, while the mode-coupling portion of the memory 
kernel describes the decay of the memory kernel at in- 
termediate and long times, and is approximated using a 
quantum mode-coupling approach. 

The fast decaying binary term determined from a 
short-time expansion of the exact Kubo transform of the 
memory function is given by 

K^f{q,t) ^ K-{q,0)e^^{-{t/T-{q)f}. (9) 
The lifetime in Eq. is approximated 



'(q) = [3if«(g,0)/4]-i/2. 



(10) 



In the above equations K'^{q, 0) is the zero time moments 
of the memory kernel given by 



^"(9,0) = 



(11) 



where uj'^"{q) is given in Eq. 
the exact expression for the 



||7J). Note that 
lifetime involves 
higher order terms such as wjj(<7) and is given 
by r«(g) = [-^'^(g,0)/2i4r«(g,0)]"i/2, where 
K-{q,0) = ~u;Uq)/u^Uq) + (oji{q) / ul{q)f . Al- 
though the calculation of (^^(9) possible it involves 
higher order derivative of the interaction potential 
and thus becomes a tedious task for the path integral 
Monte Carlo technique. Interestingly, we find that 
the approximation to the lifetime given by Eq. H10|l is 
accurate to within 5 percent for classical fluids. 3-- Thus, 
we use this simpler approximate lifetime also for the 
quantum systems studied here. 

The slow decaying mode-coupling portion of the mem- 
ory kernel, K^(q,t), must be obtained from a quan- 
tum mode-coupling approach. The basic idea behind 
this approach is that the random force projected cor- 
relation function, which determines the memory kernel 
for the intermediate scattering function (cf. Eq. ||SJ|), de- 
cays at intermediate and long times predominantly into 



modes which are associated with quasi-conserved dynam- 
ical variables. It is reasonable to assume that the decay 
of the memory kernel at long times will be governed by 
those modes that have the longest relaxation time. The 
slow decay is basically attributed to couplings between 
wavevector-dependent density modes of the form 



k,q-k — Pk Pq-k, 



(12) 



where translational invariance of the system implies that 
the only combination of densities whose inner product 
with a dynamical variable of wavevector — q is nonzero 
is pk and pq_k- 

After some tedious algebra the slow mode-coupling 
portion of the memory kernel can be approximated byi^ 



Kniq.t) 



2TT)^nJ'^{q) 
x[F'^{k,t)F^{\q^-k\,t) 
FC{k,t)FC{\<l-k\,t)] 



dk\V'^{q,k)\ 



(13) 



where n is the number density. In the above equation 
the binary term of the Kubo transform of the intermedi- 
ate scattering function, Fl^{q, t), is obtained from a short 
time expansion of F'^ (q, t) similar to the exact expansion 
used for the binary term of K'^{q, t), and is given by 



FCiq,t)^S'^{q)exp{--iuiiq)t 



1 



(14) 



The subtraction of the product of terms in Eq. (|13() in- 
volving Fl^(q, t) is done to prevent over-counting the total 
memory kernel at short times, namely, to ensures that the 
even time moments of the total memory kernel are exact 
to forth order in time. 

The vertex in Eq. H13|) is formally given by 



^''(q,k) 
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' (|i^PkPq-k) 

S'«(fc)5«(|q-k|) 

J"(g) (M/^kPq-k) 

S'^iq) S'«(fc)S"=(|q- 



(15) 



involving a double Kubo transform. "^^ In the applications 
reported below we have calculated the vertex using the 
following approximations for the three point correlation 
functions. For {PqPkPq^w) ^^'^ (Kubo) convolution ap- 
proximation has been developedfS leading to: 



1 

N 



(P^PkPq_k) 



5(fc)5'^(|q-k|)5«(g), (16) 



while for (gjjq, Pk/'q-k) have used the fact that the 
Kubo transform J'^{q) = jfijqijqi^)) can be approx- 
imated by kBT/m, m being the mass of the particle, 
within an error that is less then one percent for the rel- 
evant q values studied in this work. Based on this fact, 
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we approximate {fj^, /o£/5q_k) by: 



where t,T > 0, and 



N 



^Jqj PkPq-k 



(q-k^'^dq-ki: 



mq 



(17) 



The combination of these approximations lead to a sim- 
pUfied vertex given by 



y«(q,k) 



iksT / q-k (q — k) • q 



2mq yS'^ik) S'«(lq-k|) 



q 



Sjk) 
S^{k) 



(18) 



where S{k) is the quantum mechanical static structure 
factor (non Kubo transformed version). At high values 
of k the vertex should decay to zero, while the approxi- 
mation fails to do so. Hence, in the applications discuss 
below we employ a cutoff to overcome the shortcoming of 
our approximation. Below we discuss the choice of cutoff 
for the two model systems studied. 

To obtain the Kubo transform of the intermediate scat- 
tering function, one requires as input the frequency factor 
Lo^ijl), and the memory kernel K'^{q,t). Since the mem- 
ory kernel depends on F'^{q,t), the equation of motion 
for the intermediate scattering function (Eq. jSJ) must 
be solved self-consistently. The time-independent terms 
in the memory kernel, and the frequency factor can be 
obtained from static equilibrium input using a suitable 
path integral Monte Carlo schemeiiSiS 



III. ANALYTIC CONTINUATION OF THE 
INTERMEDIATE SCATTERING FUNCTION 

An alternative approach to the quantum mode- 
coupling theory is based on the numerical maximum en- 
tropy analytic continuation approach, which has recently 
been used by Boninsegni and Ceperley to study density 
fluctuations in liquid helium. In this section we pro- 
vide a short outline of maximum entropy analytic con- 
tinuation approach applicable to study collective density 
fluctuations in quantum liquids. 

The analytic continuation of the intermediate scatter- 
ing function is based on the Fourier relation between 
S{q,uj) and F{q,t): 



F{q,t) 



1 

2^ 



duje-"^'S{q,uj). 



(19) 



The dynamic structure factor is thus analogous to the 
spectral density used in the analytic continuation of spec- 
tral line shapesi^ and to the frequency dependent rate 
constant or diffusion constant used in analytic continu- 
ation of rates. By performing the replacement t 
—IT, and using the detailed balance relation S{q, —lo) = 
S{q,uj) we obtain 



zvr 



duj 



-f3)u 



Siq,Lo), (20) 



F(g r) = - — Tr 



(21) 



The reason for introducing the imaginary time interme- 
diate scattering function, F{q,T), is that, unlike its real 
time counterpart, it is straightforward to obtain it us- 
ing an appropriate path-integral Monte Carlo simulation 
technique^iii^ However, in order to obtain the dynamic 
structure factor and the real time intermediate scatter- 
ing function one has to invert the integral in Eq. H20() . 
Due to the singular nature of the integration kernel the 
inversion of Eq. H20|l is an ill-posed problem. As a con- 
sequence, a direct approach to the inversion would lead 
to an uncontrollable amplification of the statistical noise 
in the data for F{q,T), resulting in an infinite number 
of solutions that satisfy Eq. (|20|l . Clearly, in this case, 
little can be said about the real time dynamics and the 
corresponding dynamic structure factor. 

In recent years, Bayesian ideas have been used to deal 
with the ill-posed nature of continuing the noisy imagi- 
nary time Monte Carlo data to real timei^^ One of the 
most widely used approaches is the maximum entropy 
method. ^■^■^ The method requires only that the transfor- 
mation which relates the data and the solution be known. 
Furthermore, maximum entropy allows the inclusion of 
prior knowledge about the solution in a logically con- 
sistent fashion. As such, the method is well-suited for 
solving ill-posed mathematical problems. 

In the language of the maximum entropy method 
F(g, r) is the data, K{t, lu) = e"'^^ + e^^-'^)'^ is the sin- 
gular kernel, and S{q, oj) is the solution, also referred to 
as the map. Maximum entropy principles provide a way 
to choose the most probable solution which is consistent 
with the data through the methods of Bayesian inference. 
Typically, the data is known only at a discrete set of 
points {tj }, and thus we search for the solution at a dis- 
crete set of points {loi}. The maximum entropy method 
selects the solution which maximizes the posterior prob- 
ability, or the probability of the solution S(q,uji) given a 
data set F{q, Tj). The posterior probability is given hy^ 

V{S(q, Lo)\F{q, t)) cx exp(a,5, - x^/2). (22) 

Here is the standard mean squared deviation from the 
data, and Sq is the information entropy. 

Obtaining the maximum entropy solution then involves 
finding a map S{q,u>) which maximizes the posterior 
probability and is therefore a maximization problem in 
M variables, where M is the number of points {w;} at 
which the solution is evaluated. The solution obtained 
in this way is still conditional on the arbitrary parameter 
aq, which can be interpreted as a regularization parame- 
ter. In this work, we use a flat default map that satisfies 
a known sum rule, such as the integral over S{q,uj), and 
aq is selected according to the L-curve method^^^ In this 
case we regard aq as a regularization parameter control- 
ling the degree of smoothness of the solution, and entropy 
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as the regularizing function. The value of a, is selected 
by constructing a plot of \og[—Sq{S{q,ui))] vs. logx^- 
This curve has a characteristic L-shape, and the corner 
of the L, or the point of maximum curvature, corresponds 
to the value of aq which is the best compromise between 
fitting the data and obtaining a smooth solution. 



IV. RESULTS 

Although it is known that liquid ort/io-deuterium and 
liquid para-hydrogen may be treated as Boltzmann parti- 
cles without the complexity of numerically treating par- 
ticle statisticSjtii they still exhibits some of the hall- 
marks of a highly quantum liquid. In fact, recent theo- 
reticali4ii''^ii6.^7ii'^i^ii'^3,40,46,47,48 and experimental stud- 

ie3i2iSfliSiiSSi2& show that these dense liquids are charac- 
terized by quantum dynamical susceptibilities which are 
not reproducible using classical theories. Thus, these liq- 
uids are ideal to assess the accuracy of methods developed 
for quantum liquids such as the self consistent quantum 
mode-coupling theory and the maximum entropy ana- 
lytic continuation approach. 

Here we report on a direct comparison between these 
methods for collective density fluctuations in liquid 
ori/io-deuterium and liquid paro-hydrogen. A compari- 
son between the two approaches has been made for self- 
transport in liquid para-hydrogen^!'? where good agree- 
ment has been observed for the real time velocity au- 
tocorrelation function. However, the present compari- 
son is more challenging since the experimental dynamic 
structure factor for these liquids is characterized by 
more than a single frequency peak (unlike the case of 
self-transport),^"-^'' indicating that more than a single 
timescale is involved in the relaxation of the density fluc- 
tuations. So far, in the context of our molecular hydro- 
dynamic approach, density fluctuations in these dense 
liquids has been described only within the quantum vis- 
coelastic modeli^iiii^ which fails to reproduce the low 
frequency peak in S{q,Lj) that is associated with long 
time dynamics. Therefore, one challenge for the im- 
proved self-consistent quantum hydrodynamic approach 
and for the analytic continuation method is to recover 
this long time behavior. 



A. Technical details 

To obtain the static input required by the quantum 
mode-coupling theory and the imaginary time intermedi- 
ate scattering function required for the analytic contin- 
uation approach we have performed PIMC simulations 
at T = 20.7K and p = 0.0254A"^ for hquid ort/io- 
deuterium^S ^nd T = UK and p = 0.0235A"^ for hq- 
uid para-hydrogen.« The PIMC simulations were done 
using the NVT ensemble with 256 particles interacting 
via the Silvera-Goldman potential )^°i^^ where the entire 



molecule is described as a spherical particle, so the po- 
tential depends only on the radial distance between par- 
ticles. The staging algorithm^ for Monte Carlo chain 
moves was employed to compute the numerically exact 
input. The imaginary time interval was discretized into 
P Trotter slices of size £ = (3/P with P = 20 and P = 50 
for liquid ort/io-deuterium and liquid para-hydrogen, re- 
spectively. Approximately 2 x 10^ Monte Carlo passes 
were made, each pass consisted of attempting moves in 
all atoms and all the beads that were staged. The ac- 
ceptance ratio was set to be approximately 0.25 — 0.3 for 
both liquids. 

The static input obtained from the PIMC simulations 
was then used to generate the memory kernel {K'^{q, t) = 
Kj{q, t)+K'^{q, t))) and frequency factor needed 
for the solution of the equations of motion for F'^{q,t). 
To obtain K!^^^ {q, t) we solved Eq. with a cutoff in 
|(7— A:| to overcome the divergent behavior discussed above 
(cf. Eq. H18|) '). The choice of this cutoff is straightforward 
given that the approximate vertex decays to zero at in- 
termediated values of k before the unphysical divergent 
behavior steps in. We have used \q — k\cut — 5.66A~^ 
and \q — k\cut — 4.73A~^ for ori/io-deuterium and para- 
hydrogen, respectively. Since the memory kernel depends 
on the value of the F"^ (g, t) itself the solution must be 
obtained self-consistently. The initial guess for the mem- 
ory kernel was taken to be equal to the fast binary por- 
tion. The integro-differential equations were solved using 
a fifth-order Gear predictor-corrector algorithm.^- Typ- 
ically, less than 10 iterations were required to converge 
the correlation function, with an average error smaller 
than 10~^ percent. 

The same PIMC simulation runs were used to gener- 
ate the imaginary time F'^ (q, r) and the corresponding 
covariance matrices. We then used the L-curve method 
to determine the optimal value of the regularization pa- 
rameter, aq. In all the results shown below the value of 
aq ranged between 5 to 20 depending on the value of q 
(the results were not found to be very sensitive to the 
value of aq within a reasonable range aq = 1 — 50). The 
plots of log[—Sq{S{q, u)))] vs. logXq for all values of q re- 
sult in a very sharp L-shape curves, indicating the high 
quality of the PIMC data. 



B. Liquid ort/io-deuterium 

The results for the Kubo transform of the intermediate 
scattering function for liquid ort/io-deuterium obtained 
from the QMCT and the MaxEnt method are shown in 
the left panels of Fig. ^ for several values of q. The right 
panels of Fig. ^ show the binary and mode-coupling por- 
tions of the memory kernel obtained from the QMCT for 
the same values of q. The MaxEnt results were generated 
from S{q,oj) by taking the Fourier transform of S'(g, oj) 
divided by the proper frequency factor to "Kubo" the 
real time correlation function. 

The most striking result is the discrepancy between 
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FIG. 1: Plots of the static structure factor (upper panel), the 
Kubo transform of the intermediate scattering function (left 
panels), and the Kubo transform of the quantum binary and 
mode-coupling portions of the memory kernel (right panels) 
for liquid ort/io-deuterium. The values of q indicated in the 
panels are in units of . 



the QMCT and the MaxEnt result observed at interme- 
diate and long times for all values of g < qmax, where 
Qmax ~ 2A~^ is the value of q where S{q) reaches its 
first maximum (upper panel of Fig. The agreement 
between the QMCT which is exact to order and the 
MaxEnt results at short times is expected. However, 
the MaxEnt result fails to reproduce the proper oscil- 
latory behavior observed in F'^{q,t) that gives rise to 
a peak in S{q, oj) at a finite frequency (see Fig. |21 be- 
low), signifying the existence of collective coherent exci- 
tations in liquid ort/io-deuterium. Common to both ap- 
proaches is that as q approaches qmax the decay rate of 
the Kubo transform of the intermediate scattering func- 
tion decreases giving rise to a quantum mechanical de 
Gennes narrowing of the dynamic structure factor (see 
Fig. 121 below) . In addition, we observed an increase in 
the decay rate of F'^{q,t) at high values of q, where the 
two approaches yield nearly identical results on a sub- 
picosecond timescale. The agreement between the two 
approaches at high values of q is not surprising, since the 
decay of the intermediate scattering function is rapid, 
on timescales accessible to the analytic continuation ap- 
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FIG. 2; A plot of the dynamic structure factor of liquid ortho- 
deuterium for several values of q (in units of A~^). The thick 
solid and dashed lines are QMCT and MaxEnt results , re- 
spectively. The thin solid line is the result of QVM^ and 
open circles are the experimental results!-^. 



proach. 

The results shown in Fig. ^ are the first application 
of the self-consistent quantum mode-coup ling theory de- 
scribed in section ITU and suggested in Ref. Ilfl. As clearly 
can be seen in the right panels of Fig. the contribu- 
tion of the quantum mode-coupling portion of the mem- 
ory kernel is significant for all values of q below qmax, 
while at the highest value of q shown the contribution of 
K^{q,t) is negligible. This is consistent with results ob- 
tained for classical dense fluidsj^SiSiM signifying the need 
to include the long time portion of the memory kernel. 

To quantify the discrepancy between the QMCT and 
the MaxEnt result for F'^{q, t) we have calculated S'(g, u)) 
and compared the results to the experiments of Mukher- 
jee et alm^ for several values of q (the theoretical values 
of q are slightly different from the experiments due to the 
limitations associated with the constant volume simula- 
tion approach). The results shown in Fig. |21are normal- 
ized such that J_^dujS{q,uj) = 1. The agreement be- 
tween the experimental results2^iS& and the QMCT is ex- 
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cellent. In particular the theory captures the position of 
both the low and high intensity peaks and their width for 
nearly all wavevectors shown. The fact that for limiting 
cases the QMCT somewhat underestimates the width of 
the low intensity peak maybe attributed to the broaden- 
ing associated with the instrumental resolutioni^ This is 
certainly the case at qmax where the result of the QMCT, 
which is in excellent agreement with the result obtained 
from QVM (not shown) underestimates the width of 
the single low frequency peak. The overall good agree- 
ment between the QMCT and the experiments is remark- 
able since our previous study based on the quantum vis- 
coelastic approach failed to reproduce the low intensity 
peak associated with long time dynamics (see thin solid 
line in Fig.^). Thus, the inclusion of the quantum mode- 
coupling portion to the memory kernel is important to 
properly describe long time phenomena. In order to test 
the accuracy of the approximation for the vertex given by 
Eq. (|18l) we have replaced it with its classical limit The 
results for the low frequency part of S{q,uj) (not shown) 
are in poor agreement with the experiments. In particu- 
lar the low intensity peak is absent at intermediate values 
of q, and the higher frequency peak is somewhat shifted. 

Turning to the MaxEnt results for S{q,uj), it be- 
comes obvious why MaxEnt fails to provide a quanti- 
tative description of the density fluctuations in liquid 
ori/io-deuterium. It is well known that the MaxEnt ap- 
proach fails when several timescales arise in a problem. 
This is clearly the case here where the MaxEnt approach 
predicts a single frequency peak instead of two, at a po- 
sition that is approximately the averaged position of the 
two experimental peaks. Only when the dynamics are 
characterized by a single relaxation time, like the case 
at (/max, the MaxEnt approach provides quantitative re- 
sults. On the other hand, MaxEnt provides quantitative 
results for the short time dynamics as reflected in the 
width of the finite frequency peak in S{q, lu) (see also the 
discussion in Fig. 



C. Liquid para- hydrogen 

In many ways liquid para-hydrogen is very similar to 
liquid ort/io-deuterium. While the interaction potential is 
the same for both liquids within the Born-Oppenheimer 
approximation (for zero molecular angular momentum), 
the phase diagram is somewhat different due to the lower 
mass of para-hydrogen. Therefore, we expect that the 
accuracy of the QMCT will be similar to that obtained 
for orf/io-deuterium. As will become clear below this is 
not the case. While qualitative features are captured 
by QMCT (and not by the MaxEnt analytic continua- 
tion approach), the agreement between the experiments 
and the theory is far from the quantitative agreement 
observed for liquid ori/io-deuterium. 

In Fig. 121 we show the results for the Kubo transform 
of the intermediate scattering function for liquid para- 
hydrogen obtained from the QMCT and the MaxEnt 
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FIG. 3: Plots of the static structure factor (upper panel), the 
Kubo transform of the intermediate scattering function (left 
panels) , and the Kubo transform of the quantum binary and 
mode-couphng portions of the memory kernel (right panels) 
for liquid para-hydrogen. The values of q indicated in the 
panels are in units of . 



method (left panels) for several values of q. The right 
panels of Fig. |2| show the binary and mode-coupling por- 
tions of the memory kernel obtained from the QMCT for 
the same values of q. The agreement between the QMCT 
and the MaxEnt result is good at short times. Both 
approaches capture the quantum mechanical de Gennes 
narrowing associated with the decrease in the decay rate 
of the Kubo transform of the intermediate scattering 
function as q approaches qmax- In addition both cap- 
ture the increase in the decay rate of F'^(q,t) at high 
values of q where the two yield nearly identical results. 
However, similar to the case of liquid ort/io-deuterium, at 
intermediate and long times when q < qmax the results 
deviate markedly from each other. 

The QMCT results show pronounced oscillatory be- 
havior in F'^{q,t) that gives rise to two peaks in S{q,Lu) 
at a finite frequency as can be seen in Fig. 21 Only one 
peak is observed in the MaxEnt approach. Similar peaks 
were also observed experimentally by Bermejo et alw^ 
and using the centroid molecular dynamics (CMD)^ ap- 
proach by Kinugawa.'*'' As can be seen in the figure, 
the experimental result for the higher frequency peak is 



8 




10 20 10 20 

ffl (meV) CO (meV) 



FIG. 4: A plot of the dynamic structure factor of liquid para- 
hydrogen for several values of q (in units of A~^). The solid 
and dashed lines are QMCT and MaxEnt results , respec- 
tively. Open circles are the experimental results of Bermejo 
et al^ 

slightly shifted to lower frequencies with a smaller am- 
plitude. This discrepancy between the QMCT and the 
experiments is not completely clear. One missing factor 
is the instrument response function, which we have not 
included. We would expect that this will lower the am- 
plitude of the peaks that we have calculated. It is also 
possible that, along with inaccuracy incurred by our ap- 
proach, there may be some difficulties in extracting the 
experimental responsei^ It is interesting to note that the 
discrepancy of peak location is similar to that produced 
by the CMD result. However, the use of CMD for this 
problem is unjustified because the density operator is not 
a linear function of phase space coordinates. 

Finally, we would like to note that the agreement be- 
tween the QMCT and the MaxEnt method is excellent for 
the self-transport in para-hydrogen,^" while the results 
for the dynamic structure factor plotted in Fig. ^ show 
significant deviations between the two approaches. Most 
likely the discrepancy results from multiple frequency 



peaks in S{q,u!). MaxEnt is nearly quantitative as long 
as S{q,uj) is characterized by a single peak. Namely, at 
high q values where the decay of F'^{q,t) is on a sub- 
picosecond timescale, and at q values near qmax- 



V. CONCLUSIONS 

The problem of calculating collective density fluctua- 
tions in quantum liquids has been revisited. Two tech- 
niques have been applied to study the dynamic struc- 
ture factor in liquids orf/io-deuterium and para- hydrogen: 
The self-consistent quantum mode-coupling theory and 
the numerical maximum entropy analytic continuation 
approach. We find that the results obtained using the 
QMCT for collective density fluctuations are in excel- 
lent agreement with the experiments on liquid ortho- 
deuterium for a wide range of q values. On the other 
hand the results obtained using the MaxEnt approach 
deviate from the experiments for nearly the entire rele- 
vant wavevector range. Failure of the MaxEnt result was 
attributed to the presence of two peaks in S{q,u!). Im- 
provements of the MaxEnt, such as including real time 
dynamics in the inversion of the singular integral, may 
provide more accurate results. 

The excellent agreement between our QMCT and the 
experiments for liquid orf/io-deuterium is an encourag- 
ing indicator of the accuracy of the QMCT. The fact 
that our QMCT captures the position of both the low 
and high intensity peaks in the dynamic structure factor 
and their width for nearly all wavevectors studied is an 
important result, since our previous studies in which we 
invoked a single relaxation time for the memory kernel 
of the QGLE (the QVM) failed to reproduce the low fre- 
quency peakii^ii2ii& It is the self-consistent treatment of 
the mode-coupling portion of the memory kernel that ac- 
counts for a proper description of the intermediate to long 
time dynamics is these systems, and the proper treatment 
of the vertex in the memory kernel. 

The agreement between the QMCT and the experi- 
mental results for liquid or</io-deuterium emphasize the 
need for new experiments on liquid para-hydrogen. But 
liquid para-hydrogen is not the only system that can be 
addressed using the QMCT. Future work in other direc- 
tions, including the dynamic properties of liquid helium, 
is currently underway. 
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